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Abstract. 

We carry out a theoretical study of the collective spontaneous emission 
(superradiance) from an ultrathin film comprised of three-level atoms with V- 
configuration of the operating transitions. As the thickness of the system is small 
compared to the emission wavelength inside the film, the local-field correction to the 
averaged Maxwell field is relevant. We show that the interplay between the low- 
frequency quantum coherence within the subspace of the upper doublet states and 
the local-field correction may drastically affect the branching ratio of the operating 
transitions. This effect may be used for controlling the emission process by varying 
the doublet splitting and the amount of low-frequency coherence. 
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1. Introduction 

After the pioneering work by Kocharovskaya and Khanin [Ij |2j, Harris [3J, and 
Scully, Zhu and Gavrielides the effects of the light-matter interactions in the 
presence of low-frequency quantum coherence between sublevels in the ground or 
excited states have received much attention. There appeared new phenomena such as 
amplification and lasing without inversion (AWI and LWI), electromagnetically induced 
transparency (EIT), etc. (see references |3 |S| for review as well as the topical issue 
of Quantum Optics [7j). Optical bistability in V-type atoms has been considered in 
reference [8j. 

Recently, it was shown that in a close system of A- type atoms, the low-frequency 
coherence may give rise to superradiance without inversion (SRWI) P fTIH fTH fT ^ fT!? | (Tij . 
In reference [15J, the SRWI of an open (i.e. in the presence of a driven field) system 
of V-type atoms was discussed. We might stress that the standard superradiance (SR) 
effect, predicted by Dicke for a collection of two-level atoms jTHj, requires an initial 
inversion of level populations (see the review by Gross and Haroche [17j and the book 
by Benedict et al. ^Sj for futher details). 

The goal of the present paper is to analyze theoretically the features of the SR of 
a close system of V-type atoms. This problem was already discussed in the eighties 
by Crubellier et al jEl I2D| and by Molander and Stroud j2] within the context of 
subradiance — collective spontaneous emission from weakly radiating collective states — 
This effect has been observed experimentally by Pavolini et al. [22j in a low- 
density gas of galhum atoms (the density of galhum atoms in the interaction region 
was about 10-^^ atoms/cm^). Here, we consider a dense system comprised of V-type 
radiators and take into account the local-field correction (LFC) to the averaged Maxwell 
field. It should be stressed that this correction was neglected in previous theoretical 
studies 113 1201 and this approach is then the novelty of the paper. As is well 
known from the studies of high-density two-level assembhes |231 Ell 123 EZl EHl ESI 
Em EH E21 E21 EH ESI ESI , the LFC affects significantly the resonant optical response 
of such systems, especially the response of ultrathin films |3Zl EHl ES EDI- LFC in a 
V-type medium in the non-cooperative regime has been analyzed in reference [41j. In 
particular, the effects of the LFC on the SR of a A-type dense system have been discussed 
in references [121 [El El 132] • We show that in the case of a dense V-type system, the 
interplay between the LFC, doublet splitting and low- frequency coherence dramatically 
inffuences the branching ratio of the operating transitions. This efect can be used for 
controlhng the SR by varying the doublet splitting and the degree of coherence between 
the doublet states. 

The outline of the paper is as follows. In section El we present the model we will 
be deahng with. We derive the truncated equations for the density matrix elements 
and electric field taking into account the LFC to the Maxwell field, within the rotating 
wave approximation (RWA) . Section El is devoted to the particular case of degenerate 
doublet in the upper state, allowing for an analytical solution of the problem. Results of 
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numerical simulations are described in section^ where we provide as well interpretations 
of the peculiarities found numerically. We conclude the paper in section El with a brief 
summary of the results and how these can be related to actual measurements to infer 
the main characteristics of the LFC. 



2. Model and truncated equations 

We consider a system of three-level atoms with a doublet in the upper state (the so-called 
^-configuration, see figure QJ, forming an ultrathin film of thickness L smaller than the 
emission wavelength inside the film. This constraints is usually met in experiments and 
simplify the mathematical description of the model by taking away the propagation 
effect along the film normal. All the vectors, the SR electric field and transition dipole 
moments of the operating transitions, d2i and dai, are considered to be parallel to each 
other as well as to the film plane (the transition between the levels of the doublet is not 
included). Under these assumptions, the vector nature of the above quantities is not 
important, and the problem is reduced to its scalar form. It is also assumed, without 
loss of generality, that the transition dipole moments are real and positive: = di2 > 
and dsi = dis > 0. 

We will consider the semiclassical description of the optical dynamics of the system, 
treating the atom dynamics by means of the density matrix pap (o^,/^ = 1,2,3), i.e., 
quantum mechanically, while the field evolution is accounted for the Maxwell equation, 
i.e., classically. The effect of quantum coherence on the spontaneous emission of a single 
atom is not taken into account since is much weaker than the collective effect. 

Under the limitations adopted above, the set of equations for the joint system of 
the density matrix- Maxwell equations reads 

.dsiS' . .d2iS' , . 

psi = - luj^ip^i - i^—[p33 - Pn) - P32 , (Ij 

p21 = - lUJ21p21 - i^—[P22 - Pn) - i^—p23 , (^j 

.d2i^' , .dsiS' , . 

P32 = - ^^32P32 - ^PSI + Pl2 , (oj 

P33 = i^^{Pl3 - P3l) , (4) 
P22 = i^^^{pi2 - P2l) , (5) 

.d2iS' . .d3i8' , . . . 

Pn = - ^^-{pl2 - P2i) - ^^-{pl3 - p3i) ■ (6) 

Here, the dot denotes time derivative and S' stands for the acting field 

^' = V + —V , 7 

c 6 

where L and c denote the film thickness and the speed of hght, respectively, and 
V = No{d3ip3i + d2iP2i + c.c.) is the electric polarization of the unit volume, with 
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Nq being the atom number concentration. The first term in equation (jTj) represents the 
Maxwehian emission field, while the second term is the LFC. 

In order to further specify the model we are dealing with, we first stress that it 
is applicable to the description of SR of a thin dielectric crystalline film rather than 
a dense gas system. In the latter case, the pressure broadening terms turn out to be 
important. Having the same nature and order of magnitude as the LFC j2Z|, these terms 
have to be added to the equations for the off-diagonal density matrix elements. In solid 
crystalline media, the levels are usually broadened due to crystal imperfections as well 
as coupling to phonons. Under specific conditions, the width of the levels is smaller 
than the LFC (in frequency units). The existence of Frenkel exciton states in dielectric 
sohds, which are due to the interatomic dipolar couphng 13211111 or, in other words, due 
to the LFC, represents an unambiguous confirmation of this fact. Because of that, we 
do not take into account either relaxation of populations or dephasing of the electric 
polarization of a single atom, assuming that the SR process is faster (the estimates of 
the corresponding constants are presented in sectional). 

We seek a solution of equations dl])-® in the form: S' = exp{—iujct) + c.c, 
Psi = R31 exp{-iujct), P21 = R21 exp{-iujct), where ujc = {0J31 + a;2i)/2; and Rsi, R21 
are the complex slowly varying (in the scale 27v /ujc) amplitudes of the field and of the 
off-diagonal density matrix elements, respectively. Hereafter the latters will be referred 
to as optical coherences. Within the RWA, the equations for the amplitudes read 

^31 = - ^"|^^31 + (^-^ - i^L^ [/^3l(P33 - Pll) 

+ /i2lP32](/i21^21 + /i31^3l) , (8) 
^21 = ^^^21 + ( ) [/i2l(p22 - Pll) 

+ P3lP23](P21^21 + PsiRsi) , (9) 

Vl_ 

( 1 



P32 = -^^32P32- (^— + ^Al^ P21^3l(P21^21 +P31^3l) 

+ - ^Al^ P3li^;i(p2li^21 + P3li^3l) , (10) 

P33 = P31 + ^^^) (P21^21 + P31^3l)^31 + C.C. , (11) 

P22 = P21 {^-^ + ^^^) (P21^21 + P31^3l)^21 + C.C. , (12) 



2 

Pll = — |p2li^21 +P31^3l| • (13) 

Here we have defined /isi = d-i^jd and /i2i = d^ild, where d = {d\i + t^ii)/^; 
Al = And'^No/Sh; Tr^ = 2nkcLd'^No/h, kc = uJc/c. When deriving equations 0- 
([T^ . we exploited the fact that the equation for the slowly- varying field amplitude E^ 
can be cast in the form 

^ = + Az.) (P2ii^2i + P3ii^3i) , (14) 
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and we introduced this expression directly in the density matrix equation. The quantities 
and Al represent the magnitudes of the SR field and of the LFC (in frequency 
units), respectively [39 1. Recall that A^^ > r^^ since the relationship kcL < 1 holds for 
an ultrathin film. 

It is to be noted that equations (|H|)- (fT^ have the following integrals of motion: 

Pll +P22 + P33 = 1 , (15) 
Pu + P22 + P33 + 2(|p32r + |i^3ir + |i^2ir) = const , (16) 

where the first equation establishes the normalization condition for the total level 
population, while the second one cannot be interpreted in a simple way. 

To complete the mathematical formalism we should specify the initial conditions 
for equations (|H|)-C3|). We assume that the doublet states are initially populated, i.e., 
there exist nonzero ^33(0) and ^22(0). We also allow an initial low- frequency coherence 
^32(0). In order to trigger the emission process, we set a fixed (not-fluctuating) value for 
the initial electric polarization in the operating channels, -^31(0) = -^21(0) = Ro- This 
corresponds to triggering the SR by an ultrashort external pulse of a small area, with a 
duration Tp < min{27r/ct;32, r^^} UHl SSI SZl • 

Equations (|H|)-C3|) are written within the original basis of states |2) and |3). 
From physical reasons, especially in the case of a degenerated doublet (see below), 
another set of states turns out to be very useful: |+) = (l/\/2)(/i2i|2) + /i3i|3)) 
and |— ) = (l/\/2)(/i2i|3) — /i3i|2)). The convenience of this set is clear from the fact 
that only the superposition |+) is coupled to the ground state |1) (it will be referred 
to as bright state hereafter), while the remainder one is decoupled (dark state). The 
dipole moments of the transitions |1) |+) and |1) |— ) are = \/2d and 

) = 0. In this regards, the dark channel does not contribute to the SR. 

Within the new basis, |+), |— ), the density matrix elements can be expressed 
as follows 

^+1 = ^(/i21^21 + /i31^3l) , (17) 
P++ = ^(/^2lP22 + /^3lP33 + 2/i2l/^3l3^P32) , (18) 
^-1 = ^(/i21^31 - /^31^2l) , (19) 
P— = ^(/^2lP33 + pIiP22 - 2/i2l/i3l3^P32) , (20) 
P+- = ^[/^2l/i3l(p33 - P22) + /i2lP23 " pllP32] , (21) 

where now and p stand for populations of the bright and dark states, respectively; 

p+_ represents the low-frequency coherence, while R+i and R-i describe the coherence 
of the bright and dark channels, respectively. We stress that determines the field, 
as seen from equation The equations for these matrix elements read 

.^32 



= -i— 
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p+- 



.^32 

^^P2lPsi[P+- 

— \^+l\ 1 
Tr 



P-+ 



|^+l| 5 

rR 
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(22) 
(23) 
(24) 



.^32 



4 



/ 1 



{pll - Pli)R-i + 2fi2iP3iR+i +2 iAl) R+iP-{2^) 



\rR 



= i- 



.^32 



{pll - pIi)P+- + P2lP3l{P++ - P—) 
1 

Tr 



p— = 



.^32 / x 

1 — P21PSI(P-+-P+-) ■ 



(26) 
(27) 



As can be seen from equations 
the dark one (|— ) 



, the bright channel (|+) |1)) is coupled to 



|1)) through cc^32-terms, and thus, at ujs2 = 0, the former turns 
out to be independent of the latter (see section El for more details) . At the same time, 
the behavior of the dark channel is driven by the the bright one even in the presence of 
degeneracy of the doublet states: the field terms proportional to r^^ and Al play this 
role. 



3. Degenerated doublet {ujs2 = 0) 

We first analyze the degenerated case. Then equations (|22|)-(|2Z|), describing the bright 
channel and which we are interesting in, reduce to 
d 4 

-t;\R+i\ = — ^ \R+i\ 7 (28) 
at tr 

Z =-^\R^^W (29) 
Tr 

= -^AlZ , (30) 

where the new variables read Z = {p++ — Pii)/2, and (/) is the phase of R+i. We stress 
that equations (j2HI)-(j3D|) are similar to the SR equations of an ultrathin film of two-level 
atoms ^Hl SHI ^1 EOI, but replacing tr by tr/2 and A^, by 2A/,, which is reasonable 
due to the presence of two emission channels. 

Several qualitative conclusions about the system behavior can be drawn from the 
direct analysis of equations (|2H|)-(|3n|)- First of all, the derivative d\R+i\/dt is positive 
and thus \R+i\ will give rise to the SR if Z{0) > or, in other words, if there is an 
initial population inversion between the bright and ground states: 

P++(0) = \ [pIiP22{^) + pllP^m + 2/i2l/i3l3f^P32(0)] > pii(O) . (31) 

From here it can be shown that, in order to meet this inequality, the total initial 
population of the doublet, ^22(0) + ^33(0), must by larger than the population in the 
ground state, pii(O). In other words, in contrast to the case of a close A-system 
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where the SRWI can be observed [9J, the SR in a close V-iype system requires the 
population inversion between actual levels. For the sake of simplicity, let us set 
^22(0) = ^33(0) = ^32(0) = ^ = (1/2) [1 — Pii(O)]- It corresponds to the excitation of the 
bright state |+) = (l/\/2)(/i2i|2) +/i3i|3)) with amplitude A. Then, the inequality 
takes the form A{1 + /i2i/^3i) > Pii(O). Bearing in mind that /i2i/^3i < 1, we get 
2 A > pii(O), i.e., the total population of the doublet must to be indeed larger than that 
in the ground state. 

With the substitutions Z = i^cosO, and = 5sinG, where B ^ ^(0), 

equations (|2H|)-(|3n|) can be solved analytically, and the solution is 

't-tD\ 



^(O)sech 1^-^^ 



Z = -Z(0)tanh 



^+i| = 



(32) 
(33) 



= - 4Az. r Z{T)dT = 4Z{0)Alt'j, 
Jo 



In cosh 



t-tD 



R 



— In cosh f 



RJ J 



tD = 



2Z(0) 



(34) 
(35) 



1^+1 (0)1 

As is seen, the SR pulse is characterized by a delay time to and a duration = 
ri?/4Z(0). The unique effect of the LFC on the SR from a degenerated IZ-system is the 
SR phase modulation which changes the SR frequency 

't-to^ 



n{t) = (t) = 4Z(0) Az^tanh 



(36) 



'R 



from — 4Z(0)A/, to 4Z(0)A/,. For this reason, it is similar to what is well known for 
two-level dense systems ^Sl El El ISOj , namely within the mean-field approximation we 
are in fact dealing with, the LFC does not affect the SR kinetics, but determines the 
width of the SR spectrum. 



4. Nondegenerated doublet 

We show below that the scenario of the SR from a nondegenerated IZ-system changes 
dramatically in the presence of the LFC. It is to be noticed that, concerning the SR 
from a A-system, this fact has been already mentioned in references. jTT| fT^ fT!?| fTlj . 
In order to investigate systematically the peculiarities of the SR in the case of a V- 
system, we perform the numerical solution of equations (|H|)-(II3|)- In all our calculations, 
the dipole moments of the operating transitions |3) |1) and |2) |1) are equal to 
each other, thus implying that /i2i = /i3i = 1- The initial values of amplitudes of the 
high-frequency coherences are set to -^31(6) = -^32(0) = 10~^. Time is expressed in units 
of tr. The other initial magnitudes, such that the level populations pii(O), ^22(0) and 
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/)33(0), the low-frequency coherence ^32(0), the doublet splitting ujs2 and the LFC A^, 
will be regarded as variable parameters. 

4.1. Effects of the low-frequency coherence neglecting the LFC 

We first analyze the SR of a nondegenerated y- system setting, as a first step, the LFC 
to zero (see also the discussions in references [El EDI E])- In spite of the fact that 
this assumption might be unphysical for an ultrathin film (recall that Al > tr), the 
analysis of this ideal case will help us in understanding more comphcated situations 
with Al ^ 0. Nondegeneracy means that the magnitude of the splitting i032 is larger 
than the spectrum width of the SR in the presence of degeneracy. The latter can be 
estimated on the basis of equations as r^~^ = 4Z(0)r^^ for Al = 0. Values 

of iUs2 about several units of r^^ suffice to model the outlined condition as 4Z(0) < 2. 

Figure shows the kinetics of the SR pulse and the level populations calculated 
for = 5r^^ with the following initial conditions: all the population is in the doublet 
states, /)22(0) = ^33(0) = 0.5, and, additionally, there exists a low- frequency coherence, 
^32(0) = 0.5. Within the subspace of states |+) and |— ), this corresponds to the 
excitation of only the pure bright state |+) (p++(0) = 1), while all other populations 

and coherences are equal to zero p (0) = Pii(O) = p+_(0) = [see equations (fT7|) - 

d^]. As is seen from figure 13 the SR pulse deactivates completely the state |+). All 
the population is finally transferred to the ground state as it takes place in the case of 
two-level SR under the condition of total inversion (see, for instance, references |T7( ITHj). 
The modulation of the kinetics with frequency a;32 is explained by the fact that at 
nonzero splitting, the bright state |+) is not a stationary state: it periodically (with 
frequency CJ32, i.e. rapidly in the scale of the SR) exchanges population with the dark 
state |— ). Indeed, keeping in equations (^Hj) . (|2S|) and (^Tj) only the terms with the 
density matrix elements within the subspace of states |+) and |— ) and introducing the 
notations z = p++ — p and y = i{p+- — we obtain 



These equations describe harmonic oscillations, with frequency a;32, of a vector {y, z) in 
the YZ plane. Values z = 1,-1 correspond to the total population of bright and dark 
states, respectively. The magnitude y = 2^ya+_ refiects the low-frequency coherence. 
Note that z^ -\- y'^ = const. For the initial conditions we are dealing with (^(0) = 1 and 
2/(0) = 0), approximately one half of the period of these oscillations the system remains 
in the bright state while the other half does not. This also explains why the delay time 
of the SR in the present case (t^ ~ 18t/?) is two times as large as compared to that 
time at CJ32 = 0. Indeed, using the above initial conditions in equation one obtains 
to ^ 9tr for a degenerated doublet. 

From the above discussion it is clear that changing the sign of the initial low- 
frequency coherence, i.e. setting ^32(0) = —0.5 or, in other words, p (0) = 1, will 



y = - ^32 ^ 



Z = UJ32y ■ 



(37) 
(38) 
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not affect the SR kinetics. After a half period of oscihations with frequency 0Js2, the 
previous initial conditions will be restored. 

In figure we depicted the SR kinetics calculated for the same initial conditions 
as above, except for the low- frequency coherence P32{0) was set to zero. Within the 
subspace of states |+) and |— ), they now correspond to the excitation of an incoherent 

mixture of the bright and dark states: p++(0) = p (0) = 0.5 and p+-(0) = 

[see equations (fTHj) . ([201) and ([21]) ]• One can notice significant changes in the main 
features of the SR pulse as compared to the previous case: both the delay time 
and the pulse duration increased by approximately a factor of two and, in addition, 
the doublet states remained equally populated after the SR pulse has been emitted: 
^22(00) = ^33(00) = 0.25. More specifically, the population of the bright state, i.e. 
only one half of the total population accumulated in the upper states, is transferred 
to the ground state during the SR. During the subsequent half remains trapped in the 
dark state. The increase of the delay time and the duration of the SR pulse by a 
factor of two is simply explained by the fact that, in the present case, the initial value 
Z{0) = p++(0) — Pii(O) = 0.5 is twice as small as compared to the previous situation. 
Recall that both to and tr are inversely proportional to Z{0). 

4.2. Effects of the LFC (Al 0) 

We turn now to studying the LFC effects on the SR kinetics. Therefore, in what 
follows the magnitude of LFC, A^, is regarded as a variable parameter while the doublet 
splitting, ujs2j is set to a fixed value. More specifically, we present the results of numerical 
calculations for ujs2 = 5r^^, meaning that the splitting is larger than the full width of 
the SR spectrum which is ~ 2r^^. For small magnitudes of iUs2 compared to r^^, the SR 
kinetics is well described by equations ([32|)-([3S|)- We also assume that the total initial 
population of the system, ^22(0) = ^33(0) = 0.5 and pii(O) = 0, with equal populations 
of the doublet states, ^22(0) = ^33(0) = 0.5. As we will show below, the output depends 
on the initial value of the low- frequency coherence, ^32(0), as well. We restrict ourselves 
to two limiting cases: ^32(0) = and ^32(0) = i-y/ ^22(0)^33(0) = 0.5. The former 
corresponds to the initial excitation of incoherent mixture of the doublet states |2) and 
1 3), while the latter implies the initial excitation of the pure state |+). The initial values 
of the high-frequency coherences are set to -^21(0) = -^31(0) = 10~^ as before. 

4.2.1. SR from an incoherent mixture ^^32(0) = 0) Figure jU shows the SR kinetics 
calculated for ^22(0) = ^33(0) = 0.5 (the total initial inversion) in the absence of the 
initial low- frequency coherence, ^32(0) = 0, varying the LFC magnitude, A^. As is seen 
from this figure, increasing A^^ affects drastically the SR kinetics. While A^ is smaller 
than some "critical" value, the scenario of the SR is similar to that described in the 
previous subsection, i.e., the transitions |2) |1) and |3) |1) evolve synchronously 
and the doublet states remain equally populated after the SR pulse has been emitted: 
^22(00) = ^33(00) = 0.25. For larger A^, the transition |2) |1) in fact does not evolve. 
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conserving almost all the initial population in the state |2), while the population of the 
state 1 3) is entirely transferred to the ground state |1). This explains the changes which 
occur in the SR kinetics on increasing A/,: disappearance of the oscillatory structure 
and shortening both the SR pulse duration and delay time by approximately a factor 
of two. Indeed, as |3) |1) is the only transition contributing to the emission, the 
problem is reduced to the two-level scheme with parameters rj^ = 4:Z{0)rji = tr and 
tr, = — r^lni^3i(0) ^ ITr^^, which are characteristic for the two-level SR. It worth 
mentioning that the suppression of the transition |2) |1) occurs at d2i = d^i, i.e., 
under the condition of equivalent coupling of the individual transitions to the field. 

The physics of such a behavior is as follows. A V-atom represents two transitions 
coupled to each other by the common field which includes the emission term (~ t^^) 
and the LFC (~ A^,) [see equation ([Ti|) ]. In the presence of population inversion, the 
amplitudes of optical oscillations, R21 and i^3i, grow in time. Their increments are 
equal to each other, in the absence of the LFC and under the condition /ii3 = /ii2 and 
^33(0) = ^22(0) = 0-5. However, they become different after appearing the LFC (see 
below). This makes the transition |3) |1) to evolve faster than |2) |1). The final 
state of the system depends on the relationship between R^i and R21 at those times 
when the SR pulse is already well developed, i.e., at t ^ to- If |^2i(^d)| ~ |^3i(^d)| 
then both transitions still evolve synchronously, while at \R2i{tD)\ ^ iRsiitn)] the 
initial population of level |3) (0.5 in our case) is transferred to the level |1) before the 
oscillations |2) |1) begin to built up. It makes the population inversion between levels 
1 2) and |1) equal to zero and thus prevents the superradiant evolution of this channel 
which explains the disappearance of the oscillations of the SR pulse. 

An analysis of the linear stage of the emission, i.e., keeping all the quantities equal 
to their initial values, except for R21 and R^i, provides a solid support to the above 
arguments. Equations linearized with respect to R21 and R^i and adapted 

to the conditions used in the numerical simulations (/iai = /i2i = 1, ^33(0) = ^22(0), 
Pii(O) = ^32(0) = 0), have the form 

}j w'] Rsi + - ^Al) WR21 , (39) 



^31 = 
^21 = 



.^32 /I . . 

-^— + ( 



2 Ktr 

^ " ^Az.) H^] R21 ^{■^- ^Al) WRsi , (40) 



\rR 
( 1 



where W = ^33(0) — pii(O) = ^22(0) — Pii(O). Solving these coupled equations is 
straightforward. Below, we write down the solution in the limit r^^. A/, <C ujs2'- 

R21 ^ Roe""'' , R31 ^ Roe""'' , (41) 

where 

W / Ai 
— 1t2H^— 

As seen from equation (li^ . the increment of Rsi is indeed larger than that for R21, i.e. 



A,, = i(±^- A,W) + ^ (1 T 21^^) . (42) 

V 2 / Tr \ UJ32/ 
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Rsi grows faster than R21. Thus, during the hnear stage of the SR 

l^3l| _ ^TT.2^L t 



= exp 4H^^^- . (43) 

|^2l| V ^32 Ti?/ 

Recah that the hnear solutions for R^i and R21 are vahd almost up to the SR pulse 
maximum (see, for instance, reference [42J). Then, applying this formula for t = t^ and 
equating the exponent to unity, one obtains an estimate for A^, 

AWHd ^ ^ 

which separate two regimes of the SR. At A^, < A^, both transitions evolve 
synchronously, while for the opposite sign of the inequality, the transition |2) |1) 
is blocked for the reasons discussed above. 

Concerning numerical data {0032 = 5r^^, W = 0.5, and to ~ 35ri?), equation 
yields A^ = cc^32 / 35 = (l/7)r^^. This estimate is in good agreement with the numerical 
data (see figure We stress that the estimated value of A^ is smaller than the 
half- width of the SR spectrum given by r^^. In other words, even if the LFC (the 
dynamical resonance frequency shift) are spectroscopically hidden due to the natural 
pulse broadening, it drastically affects the SR kinetics. 

It is to be noticed that the above behavior of the SR of a l/-system resembles the 
peculiarity of the SR for the A-arrangement of nondegenerated levels ^21 ^1 E2], 
namely in the presence of the LFC, all the population from the upper level of the A- 
system is transferred to the lower level of the doublet, while the higher doublet level 
remains unpopulated after the SR pulse has gone. 



4.2.2. SR from the pure state ^^32(0) = 0.5) Figure El shows the effects of the 
LFC on the SR kinetics obtained for ^22(0) = ^33(0) = ^32(0) = 0.5 or, in other 
words, when initially the pure bright state |+) is fully populated, p++(0) = 1, while 

Pii(O) = p (0) = p+_ = 0. As is seen from figure Et, the present case differs 

noticeably from the one where initially no low- frequency coherence is created (^32(0) = 0, 
see previous section). First of all, the LFC does not affect the SR delay time at all. 
Therefore, the hnear stage of the SR is unuseful here in predicting the changes in the 
SR kinetics, as they occur when the nonlinearity is already well developed. 

The changes concern the second half and final stage of the pulse. The pulse 
shows an oscillatory structure which now cannot be associated with cc^32-oscillations. 
The frequency of the oscillations grows upon increasing A^ and, in fact, refiects the 
magnitude of the latter. The populations of the doublet states also undergo antiphased 
(with respect to each other) oscillations at the same frequency as the pulse does. This 
indicates that the doublet states start to exchange the population when the SR pulse 
is developed. A qualitative interpretation of this effect is as follows. Recall that the 
LFC shifts the frequency of the transitions |2) |1) and |3) |1) by Ax,(p22 — Pii) 
and Al{p33 — Pii), respectively. Initially, these shifts are equal to each other. FigureEb 
shows that the transition |3) |1) start to develop first. This reduces the initial 
detuning, CJ32, between the transitions. As a result, the radiation, which is emitted via 
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the transition |3) |1), is absorbed by the transition |2) |1). It further reduces the 
detuning, stronger for larger A^. After that, the transition |2) |1) starts to emit 
while |3) ^ |1) to absorb, i.e. the transitions exchange their role. 

The initial population is mostly transferred to the ground state, as it takes place for 
the same initial conditions at Al = 0. However, a small part of the population remains 
trapped in the dark state |— ), unlike the case of A^, = 0. This is the reason why the 
SR pulse has a long tail. 

5. Summary and concluding remarks 

We studied theoretically the SR from an ultrathin film of V-type atoms taking into 
account the LFC to the average Maxwell field. We show that the interplay between the 
doublet splitting, low-frequency coherence (within the subspace of the doublet states) 
and LFC may significantly affect the scenario of SR. Several conclusion can be drawn 
from our results: 

(i) Under the condition of degeneration, the three-level problem is equivalent to 
that for a two level system with a renormalized SR time. The role of the LFC is also 
similar to that for the two-level problem and manifests itself as a phase modulation of 
the SR pulse. 

(ii) For a nondegenerated IZ-system, the LFC correction affects drastically the SR 
scenario, allowing to develop one of the transitions and blocking the other one under 
specific conditions. This effect may occur even if the LFC is small compared to the SR 
spectrum, i.e., when the LFC is spectroscopically hidden. 

(iii) The SR scenario is sensitive to the amount of low-frequency coherence (within 
the subspace of the doublet states) as well as to the magnitude of the total inversion, 
thus providing a way to control the SR regimes. 

To conclude, we discuss the conditions required to prove experimentally the 
predicted regimes of the V-type SR. First of all, one should look for a dense ensemble 
of radiators where the LFC is large compared to the line width. In reference j2Hl, 
ions in KC1:0^ crystals and bound I2 excitons at donor sites in CdS single crystals 
were considered suitable for observing the LFC effects. In relation to our model, one 
should bear in mind that in disordered ensembles of dipole radiators, like the case of 

centers and bound I2 excitons, the LFC ffuctuates. The average LFC drives the 
level shifts, while the ffuctuating part contributes to the dipole-dipole line broadening. 
It turns out that both shift and line width are of the same order of magnitude as 
it takes place in dense gas systems [27j. Because of this fact, these systems can hardly 
present the effects we are discussing. 

Thin films of some organic compounds, such as naphtalene and antracene, as well 
as materials containing unoccupied d 01 f orbitals, such as Cr203 or Mn02, might be 
promising materials for this task. As they are crystalline, the intermolecular dipole- 
dipole interaction (and, subsequently, the LFC) does not fiuctuate. Furthermore, at low 
temperatures, the optical excitations in these materials are Frenkel excitons ISl^SES]- 
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This fact implies that the intermolecular dipole-dipole interaction (the LFC, in other 
words) dominates over dephasing. As the density of the optically active units 
in crystals is generally high {Nq ~ 10^^ — lO^^cm"^), the SR time constant tr 
may be small compared to the dephasing time. Indeed, tr = TiXd {^TiYdPN^L = 
(87r/3)(A/'oAc)~"'^(Ac//v)ro, where Ac = 2Ti/kc and tq = 3h/4:(Pk^ is the spontaneous 
emission time of a single emitter. Let us take Ac = 5 x 10~^ cm and assume that 
the transitions are dipole allowed (tq ~ 10~^s), which is typically the case for organic 
crystals. Then, for a film thickness L = 0.1 Ac we estimate tr as being of the order of 
10 fs. The exciton absorption line width is typically about few hundreds cm~^. This 
gives 1 ps as an estimate for the dephasing time, that is, hundreed times longer than tr. 
On the other hand, vibronic structure of aromatic crystals seems suitable for forming 
a l/-configured system. The SR of high density Frenkel excitons was observed in single 
organic crystals of R-phycoerythrin molecules at room temperature [521 • The SR pulse 
was found to be phase-modulated, thus indicating the relevance of the LFC. Therefore, 
R-phycoerythrin single crystals are promising candidates to prove the effects predicted 
in this work. 
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Figure 1. Scheme of the energy levels and transitions in a V-type atom. 



Figure 2. Kinetics of the SR field, |£:| = d\E\rR/h, the level populations pu, P22 
and p33, and the low- frequency coherence ps2 calculated for zero LFC (A^ = 0) 
at fixed magnitude of the doublet splitting LU32 = 5r^^. The initial conditions are: 
P22(0) = P33(0) = P32(0) = 0.5 and i?2i(0) = i?3i(0) = 10-8. 



Figure 3. Same as in figure O except for ^32(0) = 0. 



Figure 4. Effects of the LFC on kinetics of the SR field, |£:| = d\E\TR/h, the level 
populations pn, P22 and P33, and the low- frequency coherence P32 calculated at a 
fixed doublet splitting u;32 = 5r^^. The values of the LFC (in units of r^"^ ) are given 
the panel a). The initial conditions are: ^22(0) = ^33(0) = 0.5 (total initial inversion), 
P32(0) = (no initial low- frequency coherence) and i?2i(0) = i?3i(0) = 10-^. 



Figure 5. Same as in figure 01 except for ^32(0) = 0.5. 
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